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In this paper, we solve the eigen solutions to some nonlinear spinor equations, and 
compute several functions reflecting their characteristics. The numerical results show 
that, the nonlinear spinor equation has only finite meaningful eigen solutions, which 
have positive discrete mass spectra and anomalous magnetic moment. The nonlinear 
potential and interactions yield different contributions to the total energy, and these 
components of the energy lead to different energy-speed relation. The magnitude of 
these components can be detected by elaborate experiments. The weird properties 
of the nonlinear spinors might be closely related with the elementary particles and 
their interactions, so some deeper investigations on them are significant. 

Keywords: nonlinear spinor field, anomalous magnetic moment, mass-energy relation 
PACS numbers: 29.90. +r, ll.lO.Ef, 31.15.aj, 12.20.Fv 



I. INTRODUCTION 

Since Dirac established relativistic quantum mechanics, many scientists such as H. Weyl, 
W. Heisenberg, have attempted to associate the elementary particles with the eigenstates of 
the nonlinear spinor equation[l, 2, 3, 4, 5, 6]. In 1951, R. Finkelsten solved some rigorous 
solutions of the nonlinear spinor equation by numerical simulation, and pointed out that the 
corresponding particles have quantized mass spectra[7, 8]. However these researches have 
not realized their authors' goal due to the mathematical difficulties in analyzing nonlinear 
spinor equation. 

In recent years, a great effort has been made along this line of research. The theoretical 
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proof about the existence of solitons was investigated in [9, 10, 11, 12, 13, 14]. The symme- 
tries and many conditional exact solutions of the nonlinear spinor equations are collected in 
[15]. The present work is a development of some previous works[16, 17, 18]. In this paper, 
we define some functions which reflect the properties of eigen solutions to the nonlinear 
spinor equations, and compute the typical values, then extract some important information 
from the data. The following are some general knowledge for the nonlinear spinors. 
Denote the Minkowski metric by rj^^ = diag(l, —1, —1, —1), Pauli matrices by 



^=(^^)= h h ■ (11) 





Define 4x4 Hermitian matrices as follows 

•■={(:;)■(::)}■ < " 

where fi G {0,1,2,3}, = ct and a'^ = 7^7''. In this paper, we adopt the Hermitian 
matrices (1.2) instead of Dirac matrices 7^ for the convenience of calculation. For Dirac's 
bispinor 0, the quadratic forms of are defined by 

a" = (l)+a^'(p, j = (j)+-f(f>, /3 = 0+/?0, (1.3) 

where the superscript '+' stands for the transposed conjugation. By the definition (1.3) we 
have a'^ = 0^7'^0 etc., where 0''" = 0+7*^ is the Dirac conjugation [19]. By transformation 
law of (f), one can easily check that a'^ is a contra-variant 4-vector, 7 a true scalar and $ a 
pseudo-scalar. One can construct some other covariant quadratic forms, but they are not 
independent on (1.3) for some Pauli-Fierz identities[16, 20, 21], such as a^a^ = 7^ + 

In general, the Lagrangian of the nonlinear bispinor (p with a vector potential A^^ and 
scalar G is given by[22, 23] 

C = <P^a^{id^ - eA^)<P - //7 + 1/(7, /5) - s^G 

-)^K{d^A,d^A^ - a^A^A^) - ]^\{d^Gd^G - b'G^), (1.4) 

where A'^ and G include the self and external potential, k = ±1 and A = ±1 are used to stand 
for the repulsive or attractive self interaction, e.g. A^^ stands for repulsive electromagnetic 
potential if (k = l,a = 0), but stands for attractive interactive potential similar to strong 
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interaction if {k = —l,a ^ 0). G stands for a scalar interactive potential like Higgs field, 
which is repulsive if A = 1, and attractive if A = — 1. However, in this paper, we take (1-4) 
as one system, and only the internal interactions are considered. 

If dfjV 7^ 0, the eigen solution might be absent, so we only consider the case V = V{'j) > 
is a concave function satisfying 

V^'(7)7 > ^(7), if (7>0). (1.5) 
The corresponding dynamical equation is given by 

a^ihidf, - eA^)(j) = (/ic + sG - 1/')70> (1-6) 
(9„9° + a2)A^ = Ked^, (1.7) 
{do,d" + b^)G = As7. (1.8) 

The Hamiltonian form of (1.6) is given by 

hidt(j) = H(j), H = c[eAo + a-p+{fic + sG-V')j]. (1.9) 

where p = —hiV + eA is the momentum operator. For (1.9), by the current conservation 
law, we have the normalizing condition 

\<j)\^£x = l. (1.10) 

Let J be the angular momentum operator 

j = ffxp+^hj, 7fc = diag(crfe,(Tfc), (1.11) 

then the eigenf unctions of J3 = —hid^ + |^73 are given by 

J3(l), = JshcP,, <P, = {u^,U2e^\ivi,iV2en'^e^^\ (1.12) 

where the index 'T' stands for transpose, J3 = j + G {0, ±1, ±2, ■ ■ ■ }. For all the 
eigenfunctions, J3 is commutative with the nonlinear Hamilton operator like the linear case, 
so under a suitable choice of coordinate system, the solutions of (1.9) take the following 
form [24, 25], 

2 

(t)j = {ui,U2e'^\ivi,iv2e'^')'^ exp{jipi ^it), (1-13) 

n 
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where u^, Vk{k = 1,2) are real functions independent of tp and t. 

The system (1.9) has many symmetries such as the global gauge invariance, the spatial 
reversal invariance etc. The above assumptions have removed the uncertainty caused by 
the symmetry, and the solutions are determined except for a signature. This procedure is 
quite important for the numerical solving and stability analysis. If ^(7) take the form of 
polynomials, the solutions are analytic functions of r, and then they can be expressed as the 
Taylor series of r. In the cases of = which are the only cases for free particles, we 
have the formal solution as follows 



00 



u, + U2i = Y1 ^-e-^*'"^ + Yl ^--e^^^^^ , (1.14) 

m=0 \ n=— m+1 / 

00 / m \ 

^1 + V2^ = Y1 Jmc-^^'"-^'^' + J2 ^™ne^(^"+^)^ , (1.15) 



m=0 



where Km, Jm are real free parameters determined by boundary conditions and the normal- 
izing condition, but {Amn, Bmn) are real numbers determined by (i^„, Jn) with n < m. The 
normalizing condition (1.10) becomes 

/•oo /"TT 

2tt r^dr sin ede{uj + uj + vf + vj) = 1. (1.16) 
Jo Jo 

II. PROPERTIES OF THE DARK NONLINEAR SPINOR 

The simplest case of (1.9) is the following dynamical equation 

hidt(p = H(j), H = c[d ■ p + {fic — wj)'y]. (2.1) 

Different from the linear case, the nonlinear spinor equation generally has continuous spectra 
if the restriction of (1.10) is absent, so the normalizing condition is a quantizing condition for 
nonlinear spinors, and the nonlinear coupling coefficient w is meaningful only if the solution 
satisfies the normalizing condition (1.10). 

The eigen solutions to (2.1) with spin js = ±| can be solved rigorously as follows 



^ = ig,0,tf cose, if sm6e^^fexp{-t^t), for (P = 1, js = i) 

= {0,g,ifsmee-^\-ifcosefexp{~i^t), for (P = 1, = -i) 

V = {fcos9Jsmee^\ig,0fexp{-t^t), for (P = -1, = |) 

= {fsmee-^\-fcose,0,zgfexp{-i^t), for (P = -1, js = 



(2.2) 
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where P = 1 corresponds to even parity, and P = — 1 corresponds to odd parity. For the 
above eigenf unctions, we have 



f oo 

7=P(/-f), 47r/ ig' + fydr = l. (2.3) 

Jo 

The radial equation of even parity satisfies 

if = - - Mg' - P)]9 - If- 

For the odd parity, we have 



(2.4) 



dr-^ he 

Af=—1. 
drJ he 



(2.5) 



The initial data of (2.4) and (2.5) satisfy /(O) = 0, ^(0) > 0. For (2.4) and (2.5), we have 
positive mass < m < ^ if and only if w > 0[17]. 
Making transformation 



yu + m h (a^ + l)/i r , ^, 
' ^0 = — y „ o = ^ ' P=— ' 2.6 

.= >'°' + \ „^->'°' + ^)/. (2.7) 
y 2a/zc y 2a/uc 

where a is equivalent to the spectrum, tq takes the unit of length. (2.4) and (2.5) can be 
rewritten in a dimensionless form. For (2.4) we have 



-u' = (a — -u^ + v'^)v, u{0) = Uq > 0, 

= (i - ^2 + v'^)u - ^v, v{0) = 0, 



(2.8) 



where prime stands for ^. For (2.5) we have 



u' = d + u'^ -v'^)v, u{0)=uo>0, 
v' = {a + u^ - v'^)u - -V, v{0) = 0, 



(2.9) 



The normalizing condition (2.3) becomes 

poo 2 2 

(a + a~'Y I {u' + v')p'dp = S'^ (2.10) 







where 5* is a dimensionless constant to be determined. 
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The computation shows that, for any given a > 1, there exists a sequence of initial data 
< u{0)i < u{0)2 < such that (2.4) and (2.5) have eigen solutions. The theoretical 
analysis proves that there are infinite eigen solutions for every a[ll]. In [17] we have shown 
three families of eigen solutions with even parity and the first family of eigen functions with 
odd parity. 

To describe the properties of the eigen solutions, we define the following dimensionless 
functions, which are continuous functions of spectrum a for the same family solutions. 

1. The dimensionless norm y{a) 



1, 



oo 



y^l\g{(a + a-'ri iu' + v')p'dp]. (2.11) 







For the same family of eigen solution, y is a continuous function of a. By (2.10), the 
normalizing condition is equivalent to the equation y = IgS. 

2. The dimensionless energy S{a) 



\ ( 1 f°° 

£ = — - I mc^ H — wc / 7^ ■ Aixr'^dr 

f^c^ V 2 io 



(2.12) 



^0 

This definition of energy is in the Nother's sense. 

3. The mean diameter of an eigen solution d{a) 

2 frm'x _ a' + lJ,^{u' + v')p'dp 
~ A J\<p\^d^x a J^{u^ + v^)p^dp' ^ • ^ 

where A = ^ is a universal Compton wave length for all solutions. 

4. The total dimensionless inner pressure P{a) 

P = [mc? — [ (/ic^7 H — WC7^) ■ Arcr'^dr] 

\ Jq 2 J 



1 /a^ - 1 /o°°(«2 _ ^2)^2^^ ^ J-(^2 _ ^2)2^2^^ 



3 Va' + 1 Ji^iu^ + v^)pHp a? + 1 j^{v? + v^)pHp 



(2.14) 



The physical meanings of y(a), £{o) and d{a) are evident. Now we examine the meanings 
of P{a). For the perfect fluid in relativity, the energy momentum tensor is given by 

T^'^ = (p^ + P)U^U'' - Pg^\ = p^- 3P (2.15) 
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For the static fluid, we have the 4-dimensional speed = {1, 0, 0, 0), and then 

1 



7^0 

J- n 



For the nonhnear spinor (2.1) in curved space-time with diagonal metric, we define the 
corresponding concepts as follows [26, 27, 28] 

T'^'^ = \'Sl{(f)+{QHd'' + Q''i^^')<i)) - Cg^"" 
For static spinor, we have 

^ ^0 



p = -(t°-t;). 



(2.16) 



(2.17) 



(2.18) 



The dimensionless form of the total inner pressure of the spinor becomes (2.14). 

The curves of the dimensionless functions defined above are shown in Fig.l and Fig. 2. In 
Fig.l, the normalizing condition y = IgS = 0.918 and y = IgS = 0.647 are derived from 
the anomalous magnetic moment (AMM) of an electron according to different definition of 
mass, as computed in the next section. A rough computation was once given in [IS]. 




1 1.5 2 

Dimensionless Spectra lg(a) 



Figure 1: The norm function y{a), dimensionless energy £ = ^ and mean diameter d{a) of a 
spinor. Only the solutions corresponding to the intersection y{a) = Ig 5" are meaningful in physics 

For an electron, we have /i=me = 9.11 x lO^^^kg, H = 1.055 x 10~^^J.s, c = 2.998x lO^m/s. 
By (2.10) and S = 8.277, we can estimate the value 



w 



=4.945 X IQ-^^S^ = 3.385 x lO-^^(Jsm^). 



(2.19) 
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Dimensionless Spectra lg(a) 

Figure 2: The total energy £{a) and inner pressure P{a) of a dark spinor 

In this case, the nonhnear spinor equation has only two vahd eigen solutions corresponding 
to a = 1.95 and a = 45.7. The norm function y{a) of all other families of eigen solutions 
have no intersection points with y = 0.918. 

The radial functions {G, F) of solutions with even parity are shown in Fig. 3, where 

The unit of the coordinate r is the universal Compton wave length So the images of 
different solutions are visually comparable in Fig. 3. More images see [17]. 

III. THE NONLINEAR SPINOR WITH ELECTROMAGNETIC INTERACTION 

The nonlinear spinor with self electromagnetic interaction was researched by a few au- 
thors. In 1966, M. Wakano has approximately analyzed the cases of Aq dominance and A 
dominance when w = 0, and reached the following conclusions [29]. In the case of A domi- 
nance, the eigen solutions or the solitons do not exist for the first order approximation. In 
the case of Aq dominance, the eigen solutions exist but all with negative energy. In fact, 
the negative mass is equivalent to change the sign of Aq, which implies to transform the 
repulsive potential of Aq into the absorbent one. M. Soler and A. F. Ranada calculated the 
eigen solutions of (1-9) by omitting A. But they neglected the normalizing condition and 
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Figure 3: The radial distribution of the nonhnear dark spinors 

did not use the true value of e. Their explanation for the results seems to be misguided by 
some inadequate ideas [30, 31]. Besides, the eigen solutions with Born-Infeld potential were 
studied in [32]. The detailed non-relativistic approximation of the many-spinors equations 
was given in [17, 33] 

In general, the coordinates r and 9 can not be separable for nonlinear spinor with vector 
potential due to the term A. However Uk and Vk can be conveniently expressed by Fourier 
series with respect to 9 as (1.14) and (1.15), and the equations of the radial functions can 
be derived via variation principle, because the eigen solutions are the critical points of the 
following energy functional 

J = 2tt r^dr / sin 9d9{(p-^H(j) + -wc-f^ - -cVA^ ■ VA'' - mc^ao) + mc^. (3.1) 
Jo Jo 2 2 

So the problem (1.9) can be changed into an ordinary differential equation system which 
can be solved by numerical computation. 
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In this paper, we only consider the eigen solutions with ^-spin and even parity, which is 
the only valid case for a free electron. In the dimensionless form, we have the magnitude for 
the fields 

- a 11 

1^1 |A|~a|^|, I/I ~ -1^1, "=T^) (3-2) 

Qj Qj -Lo I 

where a is the dimensionless spectrum. Since the high order terms are caused by the vector 
potential 1^41 ~ f Ifi'l? for adequately large a, we only keep the first order approximation for 
simplicity. Then we have 

2 

0=(^, 0, if cos B, if sin Oe'^'f exp(-z-^f), (3.3) 

where g and / are real functions of r with g{0) > 0. For large spectrum a = 49.12, the 
relative error of the approximation is less than 10~^, so the approximation is accurate enough 
to reveal the anomalous magnetic moment of a spinor with electromagnetic field. The less 
the value of a, the large the error of approximation. 
The quadratic forms of (p are given by 

ao = / + /^ 7 = /-/^ a = 2^/sin6'(-sin(^,cosv2,0). (3.4) 

Correspondingly we have 

Ao = Ao{r), A = A{r) sin 9 {— simp, cos ip,0). (3.5) 

Substituting (3.4), (3.5) into (3.1) we get the energy functional 

J = Aire r^dr U[{f' + -f)g - g' f] + (^ - m)c/ - (^ + m)cf^ - -w{g^ - f^f + 

e(/ + f )Ao - \egfA + ^^(5^ + -dr)A, - \A{dl + -d^ - ^)a\ + mc\ (3.6) 

The approximation is only caused by the vector potential A. By variation, we get a closed 
system of ordinary differential equations 

g' = -i[(^ + m)c - eAo - w{g' - f')]f - l^eAg, 

f = -i[(/i - m)c + eAo - wig' - f')]g + (|.eA - f )/, (3.7) 
A'i + lA'o = -e(/ + f), A" + U' -^A= -2egf. 
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Make transformation 



h 



9 



yU + m 
/i — m' 



1 



a 



Anh 137.035999' 



r 

P = — , 



tyro 



/i'" \ h ' h ~^ 3h 

where P is dimensionless potential, which can not be confused with the pressure 
(2.14). Substituting them into (3.7), we get the dimensionless form 



(3.8) 

(3.9) 
defined in 



u - 
v' -- 
P" 



[a-P 
n , p 



a 



n2 + i;2)u + (g - p^7, 



2 pi 

p 
p 



(3.10) 



Q" + W - 



f^W+^PWdp^ 
2 4a uv 



3 J'^{u-2+v^)p2dp^ 



In (3.10) only a is a free parameter, which acts as the spectrum similar to the dark case of 
e = 0. The normalizing condition is still (2.10). (3.10) is independent on the undetermined 
coefficient w, but it becomes a global problem. The natural boundary conditions are given 

by 



u{0) > 0, v{0) = P'iO) = Q{0) = Q'(0) = 0, 
The solutions of (P, Q) can be expressed as 



(3.11) 



P 

Q 



a 



-Aa 



[u\t)+v^{t)\ T^drdp, 



p roo 

/ u{r)v{r)dTdp. 



(3.12) 
(3.13) 



We have P > 0,(5 > for the meaningful solutions. The solution of (3.12) and (3.13) can 
be soundly solved by iterative algorithm. 

The total energy of the system in Nother's sense is given by 



E = 



J2J^9J-C]d'x 



V/ 



didj) 



POO 2 1 

= 2iT r^dr / sin 9 d9 {(p^ H(j) + -wc-f^ - -cVA^ ■ VA^), 
Jo Jo 2 2 

Substituting (3.8), (3.9) into it, we get the dimensionless form 

a2-l a j;;^[{u^-v^f-P{u^ + v^)-2Quv]p^dp 



^^4 



1 



j^{v? + v'^)p^dp 



(3.14) 



(3.15) 
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The mass of a particle is a complex classical concept, which depends on the method 
of measurement and the context of theory. Using different definition of mass, we will get 
different spectrum a and constant S. Different contribution to the energy has different 
energy-speed relation, which can be detected by elaborate experiment [2 o]. Such experi- 
ment might be a key to disclose the structure of elementary particles. We give some more 
discussions in the next section. 

In what follows, we take mg and as the classical mass for computation. To get the 
anomalous magnetic moment, we introduce an infinitesimal external magnetic field 

Sext = (0,0,5), A.M = ^B{-y,x,0) = ^Brsme{-sm^,cos^,0). (3.16) 

Adding Ag^t to (3.5) and substituting it into (3.6), we get the increment of the energy 

Q roc 

AE=\—ec gfr^dr\B = n,B, (3.17) 



3 

where is the magnetic moment of the spinor. The dimensionless form is given by 

2{a^ + l) k\J^uvp^dp\ _ eh 

3^ Jo [U +v^)p^dp Znik 



where the constant ps is the Bohr magneton. 




(3.19) 



By (3.18), we get the anomalous magnetic moment of a particle 



. _ /i. -/ijj _ 2{a^ + l)k\J^uvp^dp\ 
Pb 3a jg [w^ + v'^)p'^dp 

The empirical value of the AMM of an electron is A^f = 0.001159652. The computational 
result suggests that (3.20) might be the truth of the AMM. 

To compare with the dark spinor, we also define the dimensionless norm by (2.11). The 
normalizing condition (2.10) is equivalent toy = \g S. The dimensionless functions (£", A(/, y) 
are all continuous functions of a for the same family of solutions. Fig.4 shows how to 
determine the spectrum a by the empirical AMM. Different definition of mass leads to 
different value of a. 

In Fig.4, the trends of A^f shows that A^f is a decreasing function of a, and Ag 
as a ^ oo. By the empirical data of Agf, we can compute the following undetermined 
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Figure 4: The anomalous magnetic moment of the system (3.10) vs. the spectra a, the true value 
for an electron is Ay = 0.001159652 or lg{Ag) = -2.936 

1.3 
1.2 
1.1 



0.9 
0.8 
0.7 
0.6 
0.5 

I 

Figure 5: The dimensionless functions {£{a), Ag{a),y{a)). The constant S is determined by the 
intersection points A or B, which correspond to the empirical anomalous magnetic moment 

parameters, If taking = Sfi, we have 

a = 49.12, S = 8.277, w = 3.385 x 10~^^Jsm^ Ey = 1.088keV, Ea = 85eV. (3.21) 




1.05,3=11.35 1.69,3=49.12 



0.5 



1 1.5 2 

Dimensionless Spectr3 lg(3) 



2.5 
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If taking rrik = jJ,, we have 

a = 11.35, S = 4.434, w = 9.723 x IQ-^^Jsm^ Ey = 15.08keV, Ea = 330eV. (3.22) 

Fig. 5 shows the reahstic values of some parameters such as the total energy S, the norm 
function y{a). The constants 5 or to is determined by normalizing condition y = lgS, and 
then all other parameters can be computed. By Fig. 5, we learn that, the value of a is larger 
than that of dark spinor, namely, the electromagnetic interaction increases the rest mass m 
of a spinor. 




Figure 6: The dimensionless radial functions, {u, v) correspond to spinor fields. (P, Q) correspond 
to the dimensionless potentials. 

Since is quite small, by (3.2) we learn that, if a > 10, the electromagnetic field only 

have a little influence on the eigen solution. The numerical results also show this conclusion. 
Fig. 6 shows the comparison of the dimensionless fields when a = 49.12. 
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IV. SOME INTERESTING PROPERTIES OF THE NONLINEAR SPINOR 

Form the above results and some other computation and analysis, we find some special 
but interesting properties of the nonlinear spinor. These unusual properties might have close 
relationship with the nature of the elementary particles. 

1. By P ^ and (2.14), for V = \w^'^ we find 

/•oo Y 

mc^ / (/ic^7 + -wc'f') ■ Airr'^dr. (4.1) 
Jo 2 

More calculations show that such relation also hold for other kind nonlinear potential 
y(7) satisfying V'^y — V > 0, namely we always have |P| <^ E. An interesting problem 
is whether the error is just caused by numerical approximation and P = is a rigorous 
relation generally valid for nonlinear spinors? 

2. All dimensionless energy S{a) have a similar trend S l{a oo). For large enough 
a, we always have E fic^. 

3. Taking equation system (3.10) as a developing system with respect to p, we find 
the initial value of v{0),v'{0) etc. are determined by u{0), then the eigen solutions 
satisfying boundary condition (3.11) only depend on u{0). This is a general feature 
for all eigen solutions of any nonlinear spinor, which implies the eigen equation of the 
nonlinear spinor is over-determined. This fact might be the underlying reason of the 
Pauli principle [24, 25]. 

4. For the nonlinear spinor equation with a scalar interaction 

a'^hidf.cj) = [fic + sG- V')-f(l), (a„a° + b^)G = Xsj, (4.2) 

similar to (3.12) and (3.13), G can be expressed as 

As /"°° 
Gir) = - / e-'^^-^^dr / m^e~'^^~'^d^, (4-3) 
^ Jo Jt 

so the solution to (4.2) can be soundly solved by iteration. For the AMM defined 
by (3.20), computations show that we always have Ag ~ similar to the above cases 
with electromagnetic interaction. This result implies that, it is inadequate to describe 
the strong interaction by a scalar field. 
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5. Some rough calculations show that, the AMM of a proton might be explained by the 
following nonlinear spinor with a strongly coupling vector interaction G^, 

a^ihid^ - eA^ + sG^)<j) = (Mc - V'h<P, (4.4) 
dad^A^" = ea^", + b^)G^ = sa^". (4.5) 

If /3 = ^ ~ 1, we have the following conclusions: (I). In the dimensionless form, the 
absolute values of (IGqI) |^|) are comparable with |</)| near the center, so the first order 
approximation (3.3) is invalid. (II). The first family of the eigen solutions is absent, 
namely, the fields {g{r), f{r)) have intersections with the horizontal axis near the 
origin, so such spinor has complicated interior structure. (HI). Although the solution 
has is = |, but the solution includes components with orbital angular momentum, 
which has strong influence on the AMM, so the value of \ Ag\ is much larger than that 
caused by the weakly coupled electromagnetic interaction. 

6. The energy functional of the nonlinear spinor system (3.1) is indefinite, so the sta- 
bility of the solutions is special. There are some works on this problem[34, 35, 36]. 
However, the usual treatment for the positive definite system might be inadequate for 
the nonlinear spinors. 

7. From the derivations of some previous works, we find that we could almost reconstruct 
the physical theories based on the nonlinear spinors, so some further investigations on 
the nonlinear spinors are worthwhile. 

V. THE TEST OF MASS-ENERGY RELATION 

The Einstein's mass-energy relation E = mc^ is one of the most fundamental formulae 
in physics, but it has not been seriously tested by an elaborated experiment, and only some 
indirect evidences in nuclear reaction suggested that it holds to high precision. From the 
above calculation, we found the interaction potentials of a particle will yield detectable 
effects, which lead to different energy-speed relation, which can be used as the fingerprints 
of the interactive potentials of elementary particles. So the experiment may shed lights on 
the nature of the interaction and elementary particles. 

In what follows, we give a detailed description for the experiment. The experiment only 
involves low energy accelerator of particles and measurement of speed. In this section, («, v) 



17 



stand for the speed of a particle, which can not be confused with the fields defined in (1.12) 
or (2.7). 

In Einstein's original paper [37], he derived the kinetic energy of a particle K, 

1 1 ( v^\~^ , , 

K = mc — moc , m = mo \1 -j , (5.1) 

which implies the total energy and the speed of a particle have the following simple relation 

E = m{v)c^. (5.2) 

However, this relation is based on the linear classical mechanics, and it has not been directly 
tested by elaborated experiment. There were once some indirect evidences in the nuclear 
reaction. The most accurate one is provided by S. Rainville et a/[38], which indicates that 
the mass-energy relation E = mc^ holds to an error level less than 0.00004% in the process 
of neutron capture by nuclei of sulfur and silicon resulting in 7-radiation. As pointed out 
by E. Bakhoum[39], although the authors claimed it is a direct test, it is actually a test 
for the energy conversion AE = Amc^ at low speed of the particles. As one of the most 
fundamental relation, a direct test for the original energy-speed relation (5.2) is necessary 
and significant. 

What more important is that, for the nonlinear spinors with interactive potentials such 
as electromagnetic one A'^, The detailed calculation shows the interaction terms result in 
the fine structure of the energy-speed relation [22, 23], and the fine structure can be used as 
fingerprints of the interactions. (3.21) and (3.22) show that, the nonlinear potential yields 
energy to a magnitude of IkeV, and the electromagnetic interaction to a magnitude of lOOeV, 
which can be easily detected by elaborated experiments. 

Hereafter we take c = 1 as unit of speed. The general representation of the energy-speed 
relation is given by 

^(.).^_^H-^,„^, (5.3) 



Vi 



u 



where (Mq, Mi, M^) are all constants of mass dimension, and Mq is the total static mass of 
the particles. Mi corresponds to interactions such as electromagnetic potential, M^ corre- 
sponds to the nonlinear self-interactive potential. The detailed explanation of the parameters 
see [23]. 
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Because of the little value of (Mi, M^) and the function In \/l — v?, the nonlinear effects 
can be easily concealed behind Mq. For example, when an electron get kinetic energy 30MeV. 
The corresponding speed reaches Ui = 0.99986c, but 

M„-M.+4Af, ^ 

That is to say, (5.3) is a stiff equation of the coefficients (Mq, Mi, M^). So we have to make 
some transformation to get meaningful solution. 

We propose the following experimental project to measure the coefficients (Mq, Mi, M^) 
in the energy-speed relation. The flow chart and experimental scheme are illustrated in 
Fig. 7. The particles with unit charge are produced by the particles source, and the ones 
at given initial speed uq are selected by a homogeneous magnetic fleld. By adjusting the 
radius r, we can control the initial speed of the particles uq. The series-wound accelerator 
is constructed by a set of uniform electrodes, which can be charged with high voltage V. 
When the selected particles pass though one pair electrodes, each particle receives an energy 
increment 6E = eV, which converts into its kinetic energy. If n pair electrodes are charged, 
then we get the total kinetic energy increment for each particle 

K = E{ui) - E{uo) = n5E = neV. (5.5) 

(5.3) and (5.5) establish the connection between the speed u and nV. The flnal speed u = Ui 
of the particles can be measured by the position R of the particle counter or fllm. Then we 
can determine the constants (Mo,Mi,M^) by fltting the curve f{ui,V) = as deflned by 
(5.5). 

Now we make some simpliflcation of (5.5). At flrst, we can solve the static mass Mq at 
low energy u = mi ^ c as follows. Assume the voltage V = Vq, in this nonrelativistic case, 
we have the approximation of (5.5) as follows 

neVo = ^(Mo - 2M, + M^)iul - ul). (5.6) 



Then we get 



f 

Mo = ms + 2Mi - Mj, 



Eo = m, + 2Mi- M^ + lmsul=E{uo), (5.7) 
nis = 2neVo(ul — ul)~^ , 
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Figure 7: The flow chart and experimental scheme to test the mass-energy relation 

where is the non-relativistic static mass of the particle in classical sense. Therefore, we 
only need to determine two little coefficients (Mi, M^) at high energy. 
By (5.3), (5.5) and (5.7), we have the following relation 

{2-ul)Mi - + In y^l - ^ |= £;(wi) y^l - - m,, (by (5.3) and (5.7)) 



[neV + E{uq)] \ll-ul - | , 



{neV + -rUsuD^l- «f - - 



rriaU 



l-ui 



(by (5.5) and (5.7)) 
{2Mi- M^)Jl-ul. (5.8) 



Denoting 



X 



l-ul 



U = {neV + -msU^)c ^, 



and substituting it into (5.8), we get 

x'Mi - [x + ln(l - x)] = (1 - x)U - X^s 



(5.9) 



(5.10) 



(5.10) is a linear equation of (Mi,M^), which can be easily solved by the method of least 
squares from a sequence of measured data ([/j, Xj), i = 1, 2, ■ ■ ■ , A^. Define the A^-dimensional 
vectors (a, 6, /) by, 

a = {xi,xlr-- ,Xn)^ (5-11) 
b = -1X1 + ln(l - Xi), X2 + ln(l -X2),--- ,XN + ln(l - Xn)] , (5.12) 
/ = [(1 - Xi)Ui - Xi^s, (1 - X2)U2 - X2ms, ■■■ , (1 - Xn)Un - XnitIs]. (5.13) 
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Then the solution of the least square is given by 

Ml = ^[b - {a ■ b)b] ■ f, = ^[a^b - {a ■ b)a\ ■ f, D = b ^a^ - {a ■ bf , {5.U) 

From the above computation, for an electron we have the typical order of magnitude for 
the parameters in (5.10), 

X ~ 1, U r^m.r^ IMeV, Mi ~ lOOeV, ~ IkeV. (5.15) 

The synchrotron radiation is much less than Mi and M^, so it can be omitted. A meaningful 
test strongly depends on the precision of the measurement data {xi, Ui), which should be of 
relative errors less than 10~^. How to promote the precision of the measurement is the key 
for the success of a test. 

Some possible solutions and its implications: 

1. If M-y = and Mi = 0, which means the Einstein's mass-energy strictly holds, and 
the particles can not be described by the classical fields. 

2. If M^ = and Mi ^ 0, this kind of particles has not nonlinear self-interaction, and 
the balance of the particles should be explained by scalar and vector interactions. 

3. If M^ 7^ 0, which means the standard model of particles is incomplete, and some 
calculations in quantum field theory should be modified. 

So no matter what result the experiment provides, the implication is always important and 
fundamental. Although the standard model of particles has achieved a lot of progresses in 
explanation of the behavior of micro particles, it is essentially a phenomenology. So the test 
of the energy-speed relation might be a shortcut to disclose the secrets of the fundamental 
particles and their interactions. 

VI. DISCUSSION AND CONCLUSION 

we have solved the particle-like eigen solutions to some nonlinear spinor equations, and 
computed several functions which reflect their characteristics. The numerical results show 
that, the nonlinear spinor equations have positive discrete mass spectra and anomalous 
magnetic moment. The nonlinear potential and interactions yield different contributions 
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to the total energy of the system, which can be used as fingerprints of these terms. The 
magnitude of the components can be easily detected by elaborate experiments. The weird 
properties of the nonlinear spinors might be closely related with the elementary particles 
and their interactions, so some deeper investigations on them are significant. 
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